Random matrix model for quantum dots with interactions and the conductance peak 

spacing distribution 
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We introduce a random interaction matrix model (RIMM) 
for finite-size strongly interacting fermionic systems whose 
single-particle dynamics is chaotic. The model is applied to 
Coulomb blockade quantum dots with irregular shape to de- 
scribe the crossover of the peak spacing distribution from a 
Wigner-Dyson to a Gaussian-like distribution. The crossover 
is universal within the random matrix model and is shown 
to depend on a single parameter: a scaled fluctuation width 
of the interaction matrix elements. The crossover observed 
in the RIMM is compared with the results of an Anderson 
model with Coulomb interactions. 

PACS numbers: 73.23.Hk, 05.45+b, 73.20.Dx, 73.23.-b 



The transport properties of semiconductor quantum 
dots can be measured by connecting the dots to leads via 
point contacts jjj . When these point contacts are pinched 
off, effective barriers are formed between the dot and the 
leads, and the charge on the dot is quantized. Adding 
an electron to the dot requires a charging energy Ec to 
overcome the Coulomb repulsion with electrons already 
in the dot. This repulsion can be compensated by mod- 
ifying the gate voltage V g on the dot. For temperatures 
below Ec, a series of Coulomb blockade oscillations is ob- 
served in the linear conductance as a function of V g . For 
temperatures much smaller than the mean level spacing 
A, the conductance is dominated by resonant tunneling 
and the Coulomb blockade oscillations become a series of 
sharp peaks. 

In dots with irregular shapes, the classical single- 
electron dynamics is mostly chaotic. Quantum mechan- 
ically, chaotic systems are expected to exhibit universal 
fluctuations that are described by random matrix theory 
(RMT). The distributions || and parametric correlations 
H of the Coulomb blockade peak heights in quantum 
dots have been derived using RMT, and these predictions 
have been confirmed experimentally . 

Another quantity of recent experimental and theoreti- 
cal interest is the peak spacing statistics. The peak spac- 
ing A2 can be expressed as a second order difference of 
the ground state energy £ g ™^ of the n-electron dot as a 
function of the number of electrons: 



An = £0 l + 1 ) 4- PC"- 1 ) 
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Using the constant interaction model (which ignores 
interactions except for a classical Coulomb energy of 
n 2 Ec/2), and assuming a single-particle spectrum that 



is independent of n, A2 = E n+ i — E n + Ec, where E n 
is the ro-th single-particle energy. Within this model, 
RMT suggests a Wigner-Dyson distribution of the peak 
spacings with a width of ~ A/2. However recent ex- 
periments find a distribution that is Gaussian-like and 
has a larger width |^-^|. This observation underlines 
the limitations of the constant interaction model and the 
importance of electron-electron interactions beyond an 
average Coulomb energy. Some observed features of the 
peak spacing distribution have been reproduced using ex- 
act numerical diagonalization of small disordered dots 
(n < 10) with Coulomb interactions f^,|io||. The width of 
the distribution is found to increase monotonically with 
the gas parameter r s . Analytic RPA estimates in a disor- 
dered dot for small values of r s give peak spacing fluctu- 
ations that are larger than those of RMT but still of the 
order of A PJ. Recent Hartree-Fock calculations p2 13 
of larger disordered and chaotic dots with interactions 
(up to n ~ 50 electrons) also reproduce Gaussian-like 
peak spacing distributions. 

The above studies were carried out for particular mod- 
els of quantum dots using Coulomb as well as nearest- 
neighbor interactions, and it is not clear how generic 
the conclusions are. It is also not obvious which bare 
electron-electron interaction should be taken to represent 
the experiments because of screening generated by exter- 
nal charges. It is therefore important to find out whether 
the observed statistics of the peak spacings is generic, and 
in particular whether it can be reproduced by a modified 
random matrix model. Standard RMT does not make ex- 
plicit reference to interactions or to number of particles. 
To study generic interaction effects on the statistics, we 
need a random matrix model in which the two-body inter- 
actions are distinguished from the one-body part of the 
Hamiltonian. Recently a two-body random interaction 
model (TBRIM) introduced in nuclear physics ju] was 
used together with a diagonal random one-body Hamil- 
tonian to study thermalization in finite-size systems Jig] 
and the crossover from Poisson to Wigner-Dyson statis- 
tics in many-body systems |l6| . The model explains sta- 
tistical features observed in atomic |l7j and nuclear |T§| ] 
shell model calculations. However, the Poisson statistics 
that was used as a non-interacting limit of the model 
Jl6| , p7| is not suitable for the study of dots whose single- 
electron dynamics is chaotic. 

In this paper we introduce a random interaction ma- 
trix model (RIMM) for strongly interacting Fermi sys- 
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terns whose single-particle dynamics is chaotic. With 
this model we can study generic and universal effects as- 
sociated with the interplay of one-body chaos and two- 
body interactions. In particular, we apply the model 
to study the peak spacing statistics and find a crossover 
from a Wigner-Dyson distribution to a Gaussian-like dis- 
tribution as a function of a parameter that measures 
the fluctuations of the interaction matrix elements. The 
crossover depends on both the number of particles and 
the number of single-particle orbits but becomes uni- 
versal upon an appropriate scaling of the interaction 
strength. The crossover is demonstrated in a model of 
a small disordered dot with Coulomb interactions, and 
we show that the results can be scaled approximately 
onto those of the RIMM. 

A general Hamiltonian for spinless interacting fermions 
has the form 

H = ^2 11 i + Uijkia\a)aia,k , (2) 

ij ijkl 

where hij are the matrix elements of the one-body Hamil- 
tonian and Uij-ki — (ij\u\kl) — (ij\u\lk) are the antisym- 
metrized matrix elements of the two-body interaction. 
The states \i) = a]|0) describe a fixed basis of m single- 
particle states. We define an ensemble of random matri- 
ces H of the form (|^) , where the one-body mxm matrix 
belongs to the Gaussian ensemble of symmetry class 
(3, and the matrix elements of the two-body interaction 
are real independent Gaussian variables with zero aver- 
age and variance U 2 {\U 2 ) for the diagonal (off-diagonal) 
elements 

P{h) oc e~^ Tr h2 ; P(u) oc e" Tr a2 ^ , (3) 

Eqs. (|) and ([§ define the RIMM. The parameter a 
determines the single-particle level spacing A. In the 
non-interacting limit U — 0, the random ensemble de- 
scribes the universal statistical properties of a system 
whose single-particle dynamics is chaotic. For conserved 
time-reversal symmetry h is a GOE matrix, while for 
broken time-reversal symmetry, e.g. in presence of an 
external magnetic field, h becomes a GUE matrix. The 
random ensemble for the two-body part of the Hamilto- 
nian (0) is invariant under orthogonal transformations of 
the single-particle basis. An average interaction that is 
invariant under such transformations can be included in 
the model, but it leads to a constant charging energy shift 
in A2 and does not affect the peak spacing statistics. 

The randomness of the two-body interaction matrix 
elements is motivated by the strong fluctuations of the 
Coulomb interaction matrix elements in the basis of 
eigenstates of the chaotic single-particle Hamiltonian. 
The RIMM differs from the TBRIM in its one- 

body part, which is less relevant at the high excitation 
energies considered in the earlier studies but is of crucial 
importance near the ground state. The TBRIM has a 



Poissonian statistics in the non-interacting limit, in con- 
trast to the Wigner-Dyson statistics characterizing the 
non-interacting limit of (Q). Our random interaction ma- 
trix model is suitable for describing the generic statistical 
fluctuations in quantum dots with chaotic single-particle 
dynamics and in the presence of electron-electron inter- 
actions. The model depends on three parameters: U/A, 
the number of single-particle orbits m, and the number 
of particles n. 
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FIG. 1. Peak spacing distributions P(Aa) for the random 
matrix model (JiJ) for m = 12, n = 4 and for U/A = (solid 
circles), 0.35 (open squares), 0.7 (solid diamonds), 1.1 (open 
triangles) and 1.8 (solid triangles). The one-body part h is a 
GOE. We see a crossover from a Wigner-Dyson distribution 
for (7 = (dashed line) to Gaussian-like distributions for 
U/A > 1. The solid lines are Gaussian distributions with the 
widths of the U/A = 1.1 and 1.8 distributions, respectively. 

Next we apply the RIMM(|2|) to study the peak spac- 
ing statistics in Coulomb blockade quantum dots. Peak 
spacings are computed using ([!]); i.e. the ground-state 
energy is calculated for three consecutive numbers of par- 
ticles n — 1, n and n + 1, and statistics are collected by 
generating realizations of the ensemble H. 

Typical distributions of A 2 = (A 2 - (A 2 ))/A for 
the case of conserved time-reversal symmetry (h is a 
GOE matrix) are shown in Fig. |l| for several values of 
U/A. For the non-interacting case we obtain the Wigner- 
Dyson distribution (dashed line), but as U/A increases a 
crossover is observed to a Gaussian-like distribution. The 
distributions for U/A = 1.1 and 1.8 are compared with a 
Gaussian of the same width (solid lines) . The model (^) 
does not include spin and therefore cannot reproduce the 
expected bimodal structure of the peak spacing distribu- 
tion at weak interactions. However, numerical simula- 
tions in small disordered dots indicate that this bimodal 
structure disappears already for weak interactions ]l0| . 

The standard deviation of the spacing fluctuations 
ct(A 2 ) (in units of A) is shown in the top panel of Fig. 
H vs. U/A for n — 4 and several values of m. The sta- 
tistical errors (due to finite number of samples) are also 
estimated but are smaller than the size of the symbols. 
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At zero interaction we are close to the GOE value of 
w 0.52. <t(A2) increases slowly and then more rapidly 
above U/A ~ 0.5. At strong interactions it is approxi- 
mately linear in U/A. The top inset of Fig. || shows sim- 
ilar curves of the spacing fluctuations but for constant 
number of single-particle states to = 14 and for several 
values of n. The standard deviation curve versus U/A 
depends on both m and n. However upon the linear scal- 
ing U e g = f(m, n)U/A (f(m, n) is a scaling constant) all 
curves coalesce to a single universal curve. 




FIG. 2. Top panel: the standard deviation a(A2) of the 
peak spacing fluctuations as a function of U/A for the ran- 
dom matrix model (j^) with a GOE h. Shown are curves for 
n — 4 and m = 10 (circles), 12 (squares), 14 (diamonds) 
and 16 (triangles). Top inset: similar curves for m — 14 
and n = 4 (circles), 6 (squares) and 8 (diamonds). Bottom 
panel: The same curves as in the top panel but as a func- 
tion of t4ff = f(m,n)U/A. The curves are shown by their 
corresponding solid symbols except for the reference curve 
(m = 12, n = 4) which is shown by a solid line. Similar 
scaled curves (n = 4 and m = 10, 14) are shown for the GUE 
case (open symbols) and compared with the GUE reference 
curve (dashed line). Left inset: the scaling factors f(m,n) 
as a function of n for the GOE (solid symbols) and the GUE 
(open symbols) for m = 10, 12, 14 and 16 (from top to bot- 
tom). Bottom inset: the ratio <7goe(A2)/o"gue(A2) versus 
U e g calculated using the m = 12, n = 4 data. 

To demonstrate the scaling we first choose a 'reference' 



curve, e.g., to = 12 and n — 4, which we determine ac- 
curately using 10,000 realizations at each value of U/A. 
For other values of (m, n) we use typically ~ 1000 — 5000 
realizations and find the scaling factors /(m, n) by a least 
squares fit. The bottom panel of Fig. || shows the same 
curves of the top panel (solid symbols) in comparison 
with the reference curve (solid line), but as a function of 
the scaled parameter U c s- The curves scale almost per- 
fectly within the statistical errors. Also shown are scaled 
GUE curves (open symbols) for n = 4 in comparison with 
the GUE reference curve (dashed line). 

The scaling factor f(m,n) is shown in the left inset 
of the bottom panel of Fig. || as a function of n for 
different values of to and for both the GOE (solid symbols 
with error bars) and the GUE (large open symbols) cases. 
Within the statistical errors f(m, n) is independent of the 
symmetry class, supporting the universality of our model. 
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FIG. 3. Left panel: a {An) versus U c /V for a 4 x 5 cylin- 
drical Anderson lattice with Coulomb interactions. Shown 
are curves for $ = and a disorder strength of W = 3, 5, 7 
(solid circles, squares and diamonds, respectively), and for 
$ = ClS^o and W = 5 (open squares). Right top panel: 
<j(A2)/A versus U c /V for the same cases shown on the left 
panel. Notice the weak dependence on W . Right bottom 
panel: <r(A2)/A for the Coulomb model (symbols) versus the 
scaled Ueg = coU c /V in comparison with the reference curves 
of the GOE (solid) and GUE (dashed) random matrix model. 

The width of the spacing distribution is larger for 
the GOE case than for the GUE case at any value of 
U c ff. The bottom right inset of Fig. || is the ratio 
ctgoe(A2)/ctgue(A2) versus U c s, calculated from the 
reference curves. The ratio is ~ 1.24 ± 0.02 for the 
non-interacting case (in close agreement with the RMT 
value), and depends only weakly on the interaction in 
the crossover regime U e g < 1. This is consistent with 
recent measurements [^| which find a ratio of ~ 1.2 — 1.3 
for semiconductor quantum dots with a gas constant of 
r s ~ 1 — 2. At stronger interactions the ratio decreases. 
At large values of U, the two-body Hamiltonian domi- 
nates and one can ignore the one-body part. Since it is 
only the latter that distinguishes between the conserved 
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and broken time-reversal symmetry cases, the ratio of the 
widths approaches 1 at strong interactions. 

Next we compare the crossover observed in the RIMM 
to the results for a tight-binding Anderson model with 
cylindrical geometry, hopping parameter V = 1, and on- 
site disorder with a box distribution of width W. Elec- 
trons at different sites interact via a Coulomb interaction 
whose strength is U c = e 2 /a over one lattice constant 
a. The standard deviation of the peak spacing a {£±2) is 
shown in the left panel of Fig. || versus U c /V for a 4 x 5 
lattice, n = 4 electrons and several values of W. The 
values of W are chosen so that the RMT statistics is ap- 
proximately satisfied in the non-interacting case. In the 
absence of a magnetic field we choose W = 3,5,7. How- 
ever, in the presence of a magnetic flux, which we apply 
inside the cylinder and incorporate in the hopping matrix 
elements in the perpendicular direction ($ = 0.15$o), 
only the W = 5 case satisfies the spectral RMT statistics. 
We find that 03=0(^2) is monotonically increasing with 
W. After rescaling a by the mean level spacing A at the 
Fermi energy, the residual IF-dependence of cr(A 2 )/A is 
rather weak (top right panel of Fig. ||). The standard 
deviation curves for the Coulomb model can be mapped 
approximately on the random matrix model curve (bot- 
tom right panel of Fig. ||) by defining U c s — cqU c /V for 
some constant Co that depends on the disorder strength 
and lattice size, cq depends only weakly on W. 
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FIG. 4. Peak spacing distributions for the random matrix 
model (^) (where h is a GOE matrix) for m = 10, n = 6, 
U/A = 0.24 (open circles); m = 12, n = 4, U/A = 0.33 (open 
squares), and m = 14, n = 4, U/A = 0.42 (open diamonds). 
In all three cases U e s = 0.33. The solid triangles show the 
peak spacing distribution of the Coulomb model for n — 4 
electrons on a 4 x 5 lattice with W = 5 and U c /V = 0.75. 

The universal aspects of the crossover can also be in- 
vestigated by studying the peak spacing distributions 
themselves. In Fig. ^ we show the peak spacing dis- 
tributions for three different values of (m, n) but at con- 
stant U e g = 0.33 (corresponding to three different val- 
ues of U/A). We find that all three distributions coin- 
cide, indicating that finite size effects in the random ma- 



trix model are negligible. A corresponding distribution 
for the Coulomb model is also shown for U c /V — 0.75, 
W = 5 and n = 4 and is rather close to the random 
matrix model distributions. The deviations seen in the 
Coulomb case may be partly due to finite size effects that 
are non-universal; even at U c = we observe deviations 
from the expected Wigner-Dyson distribution. 

In conclusion, we showed that a random interaction 
matrix model that includes a one-body part belonging to 
one of the standard Gaussian ensembles and a two-body 
random interaction is suitable for studying generic in- 
teraction effects on the statistics of finite Fermi systems 
whose single-particle dynamics is chaotic. We applied the 
model to chaotic dots in the Coulomb blockade regime, 
where it describes a crossover of the peak spacing statis- 
tics from a Wigner-Dyson to a Gaussian-like distribution. 
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